Contractile stresses in cohesive cell layers on finite-thickness substrates 
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Using a minimal model of cells or cohesive cell layers as continuum active elastic media, we examine 
the effect of substrate thickness and stiffness on traction forces exerted by strongly adhering cells. 
We obtain a simple expression for the length scale controlling the spatial variation of stresses in 
terms of cell and substrate parameters that describes the crossover between the thin and thick 
substrate limits. Our model is an important step towards a unified theoretical description of the 
dependence of traction forces on cell or colony size, acto-myosin contractility, substrate depth and 
stiffness, and strength of focal adhesions, and makes experimentally testable predictions. 



Many cell functions, such as spreading, growth, differ- 
entiation and migration, are affected by the elastic and 
geometric properties of the extracellular matrix [1 . Con- 
siderable effort has been devoted to the study of cell ad- 
hesion to elastic substrates [2]. Cells adhere to a sub- 
strate via focal adhesion complexes that link the sub- 
strate to the actomyosin cytoskeleton, which in turn gen- 
erates contractile forces that deform soft substrates [3]. 
The traction forces that the cell exerts on the substrate 
are regulated by the cell itself in a complex feedback loop 
controlled by cell activity and substrate elasticity. 

Two powerful experimental techniques have been de- 
veloped to measure forces by cells on substrates: traction 
force microscopy, used to probe cell adhesion to contin- 
uous substrates [H [5J, and the imaging of cell-induced 
bending of microfabricated pillar arrays [BJ. These two 
techniques have also been recently combined [7]. These 
experiments have yielded new insight on substrate rigid- 
ity sensing and have opened up new questions on the 
physics of individual and collective cell adhesion: What 
controls the length scale that governs the penetration of 
traction forces? What is the relative role of active cel- 
lular contractility and cell-cell-interaction in controlling 
the emergent response of cell layers? In this Letter we 
describe minimal models of individual cells and adhering 
cell colonies that reproduce qualitatively several experi- 
mental findings. The traction stresses exerted by cells on 
substrates are extracted directly from measurements of 
micropillar displacements or inferred from the displace- 
ments of fiducial markers embedded in a continuum sub- 
strate. It is found that traction stresses by isolated fi- 
broblasts and epithelial cells on pillar arrays are local- 
ized near the cell edge, while contractile stresses (referred 
to below as cellular stresses) built-up inside the cellular 
material is largest near the cell center [H [8], as shown 
schematically in Fig. [T] This behavior, also observed in 
adherent cell sheets and in migrating cell colonies [5HTT] , 
is predicted by our model. Further, both substrate thick- 
ness and stiffness affect cellular and traction stresses [TJ] . 
The magnitude of the traction stress increases with sub- 
strate stiffness, saturating at large stiffness [8], and it de- 
creases sharply with substrate thickness, indicating that 
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FIG. 1. Schematic of a cell layer of lateral extent L and thick- 
ness h c « L adhering to a substrate of thickness h s . The 
build-up of contractile stress a in the cell layer is indicated 
by the color map, while the traction stresses in the substrate 
are shown as vectors (blue online). The spatial variation of 
both traction and cellular stresses in the lateral (x) direction 
are characterized by the length scale £ p , referred to as the 
penetration length. 



cell colonies on thick substrates only probe a portion of 
substrate of effective depth comparable to the lateral ex- 
tent of the cell colony [T3]. Both trends are reproduced 
by our model (Fig. [3]). 

Our model builds on recent work [131 ITS] describing the 
cell or cell layer as a contractile elastic medium, with local 
elastic response of the substrate (as appropriate for mi- 
cropillar arrays or very thin substrates). In contrast, here 
we consider substrates of finite thickness where the non- 
locality of the elastic response must be included. While 
previous studies have analyzed the deformations of finite- 
thickness substrates due to point traction forces on their 
surface |16l 117] , our work considers the inhomogeneous 
traction due to an extended contractile cell layer. A cen- 
tral result for our work is the expression for the scaling 
parameter referred to as the lateral penetration length t p 
(Fig. [T]) . This length scale characterizes the in- plane spa- 
tial variations of both adhesion-induced traction stresses 
on the substrate and cellular stresses within the cell layer 
in terms of cell and substrate elastic and geometrical 
properties. Our model also quantifies the experimentally- 
observed role of substrate thickness h s in controlling the 
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FIG. 2. Top : Schematic of a contractile cell adhering 
to a soft substrate (left) and effective spring constant fc e ff 
versus cellular strain A£, showing strain stiffening (right). 
Bottom : Cell contraction A£ (solid blue line) and trac- 
tion force Ft (red dashed line) vs substrate stiffness (left) 
for h s = 10 fim and as a function of substrate thickness 
(right) for E s = 500 Pa. Other parameters : F A = 10 nN, 
k c = lnN//im, k a = 2.5nN//im, E s — IkPa, h s — 10 ^m, 
4o = 10 /im, v a = 0.4. 



mechanical response of adhering cell layers [T2] . If h s is 
small compared to the lateral extent L of the cell sheet, 
the substrate elasticity plays a negligible role in deter- 
mining the mechanical response of the cell. This may 
explain why traction forces exerted by cell colonies with 
L > h s appear insensitive to substrate stiffness |llj . If, 
in contrast, L <^ h s , then substrate nonlocality controls 
stress build-up in the cell sheet. This crossover may be 
observable in large cell colonies on thick substrates. Fi- 
nally, the importance of long-range substrate elasticity 
has also been emphasized in recent models of cells as 
active dipoles on a soft elastic matrix, where it is cru- 
cial in controlling cell adhesion [TH1 [H] . Long-range in- 
terfacial elastic stresses coupled with gel thickness have 
also been shown to have a profound effect on focal adhe- 
sion growth [2D] and to enhance cell polarization [3TJ [52] . 
These important effects are not discussed here. 

Contractile cell on a soft substrate. To illustrate the 
importance of substrate nonlocality, we first analyze a 
single cell, modeled as a contractile spring of stiffness k c 
and rest length £ c q, adhering to a continuum substrate 
(described as an elastic continuum of Young's modulus 
E s and Poisson's ration v s ) via two focal adhesion bonds 
(linear springs of stiffness k a ) located at x\ and x 2 (Fig. [2] 
top left) [23]. This is motivated by the experimental 
observation that in adhering cells focal adhesions tend to 
be localized near the cell periphery [53] . For simplicity we 
consider a one dimensional model, where the cell lies on 
the x axis and the substrate lies in the < z < h s region 



of the xz plane. Contractile acto-myosin fibers connect 
the focal adhesions and exert active forces of magnitude 
Fa- Once the cell has fully adhered, the cell-substrate 
system is in mechanical equilibrium. Force balance at X\ 
and x 2 yields 

k a [ui - u s (xi)} = Fa - k c (m - u 2 ) , (la) 
k a [u 2 - u s (x 2 )} = -F A + k c {m - u 2 ) , (lb) 

with Ui the displacements of the contact points Xi from 
their unstretched positions x 2 — x\ — £ c q, and u s (xi) the 
displacement of the substrate's surface at Xi. All dis- 
placements are defined with respect to an initial state 
where the cell has length £ c0 . The net contraction is 
then A£ = I c q — (x 2 — xi) = u\ — u 2 . The trac- 
tion force by the cell on the substrate is localized at 
x\ and x 2l yielding a traction force density frix) = 
F T S(x - xi) - F T S(x - x 2 ), with F T = F A - k c A£. As- 
suming linear elasticity, the substrate deformation is [25] , 
u s (x) — dx'G(x— x')fT(x'), where G(x) is the elastic 
Green's function at z = h s . For a substrate of thickness 
h s we use the approximate form 26 
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derived in the Supplemental Material [57], with a the 
size of adhesion complexes, providing a short-distance 
cut-off, and Kq denotes the modified Bessel function of 
the second kind. We obtain F T (A£) = ±k cS (A£)A£, with 
fc^ 1 = k~ 1 + [G(0) - G(l c0 - A£)] the effective stiffness of 
the cell-substrate adhesions. For A£ -C Z c0 , k e g is inde- 
pendent of A£ and Ft scales linearly with A£. Stiffening 
sets in for A£ > £ c q\1 — h s (l + v s )/£ ca \ : as shown in Fig. [5] 
(top right), with a crossover controlled by the thickness 
of the substrate h s . Using Ft = Fa — k c A£, we solve 
for both A£ and Ft, shown in Fig. [2] (bottom) as func- 
tions of the substrate thickness and stiffness. For very 
thin (h s — > 0) or infinitely rigid substrates, where the 
substrate elasticity becomes local, A£ = Fa/ {k c + k a /2), 
corresponding to a spring k c in parallel with a series of 
two focal adhesions springs k a . In this limit the traction 
force saturates to Ft — k a FA/(2k c + k a ). Conversely, 
for a very soft substrate with E s — > 0, the contraction 
is maximal and given by F^/fcc, and Ft — > 0. The sub- 
strate thickness above which both cell contraction and 
traction force saturate is controlled by the cell size and 
the substrate elasticity, in qualitative agreement with ex- 
periments [T2J. 

Contractile Cell Layer. The continuum limit can be 
obtained by considering a multi-mer of N — [L/I c q] con- 
tractile elemental ''cells" , connected by springs represent- 
ing cell-cell interactions. The outcome is a set of coupled 
equations for a contractile elastic medium. For a cell 
layer of thickness h c « L (Fig. [lj, the force balance 
equation, averaged over the cell thickness, is 



Y a [u(x) - u s (x)] = h c d x a(x) , 



(3) 
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where Y a — k a /(L£ c o) describes the effective strength 
of the focal adhesions, u(x) is the displacement field 
of the cellular medium at z = h s , and a is 
the thickness-averaged cellular stress tensor, o(x) — 
V^-c fh° +Hc dz o- xx (x,z), given by a(x) = B c d x u + a a , 
with B c the longitudinal clastic modulus of the cell layer. 
The one dimensional model presented here may be rel- 
evant to wound healing assays, where the cell layer is a 
strip with y-translational invariance. Although we have 
neglected components of the cellular displacements and 
spatial variations along z, the cell elastic constants are 
those of a three-dimensional cellular medium. The active 
stress a a — Fa/{LIi c ) arises from acto-myosin contractil- 
ity [28]. The substrate deformation at the surface is 



u s (x) = h c / dx'G(x - x')& x a{x') , 



(4) 



with G(x) the elastic Green's function of a substrate of 
infinite extent in x, occupying the region < z < h s , 
evaluated at z = h s . Eqs. Q-Q can be reduced to 
integro-differential equations for the cellular stress, as 



l 2 a dl<j + <j a = a-B c Lh c dl [ dx'GQx - x'\)a{x' 

Jo 



)• (5) 



The length scale l a = \/B c h c /Y a controls spatial varia- 
tions of cellular stresses induced by the stiffness of the 
focal adhesions. It is the size of a region where the 
areal elastic energy density Y a i 2 a associated with focal 
adhesions is of order of the areal elastic energy den- 
sity B c h c of the cell layer. For a cell monolayer with 
B c = 1 kPa, h c = 0.1 /im, L = 100 /im, £ c q = 10 /im and 
k a = 2.5nN//mi [3], we get t a ~ 6.3/im, comparable to 
traction penetration length seen in experiments on stiff 
microposts [29 , 30 . The second term on the right hand 
side of Eq. (|5| describes spatial variations in the cellular 
stress due to the (generally nonlocal) coupling to the sub- 
strate. In the following we examine solutions to Eq. ([5]), 
considering various limiting cases for the substrate thick- 
ness and analyze the dependence of traction stresses on 
cell size, substrate stiffness and substrate depth. The 
equation governing stress distribution in two dimensional 
cell layers is derived in the Supplemental material |27j . 

Thin substrate. If the substrate's elastic response can 
be approximated as local, as it is the case for h s « 
L or for cells on micropillar arrays, the Green's func- 
tion is given by G(x) = 2hs( ^+ v =) 5(x). Eq. ^ can 
then be written as l 2 p d\a + a a = a, where, l p = 
y/B c h c /Y eS and = Y' 1 + 2h s (l + u s )/E s describes 
the combined action of the focal adhesions and the sub- 
strate, acting like two linear elastic components in sc- 
ries. Assuming zero external stresses at the boundary, 
i.e., cr(0) = <j(L) — 0, the internal stress profile is 
a(x) = cr a (l- cosh [{L-2x)/2£ p }/ cosh [L/2£ P ]) QUI OH 
[T5] . The traction stress T(x) = Y c su(x), is localized 
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FIG. 3. Internal stress a(x)/a a (top left), substrate displace- 
ment u s (x) (top right) and traction stress T(x)/a a (bottom) 
vs position x along the cell layer, for E s = 500 Pa (solid, 
blue), 50 Pa (dotted, green) and 10 Pa (dashed, red). The 
vertical dashed lines in the top right frame denote the cell 
layer's edges. Inset (bottom right): magnitude of contractile 
moment \V\ vs E s /B c . Inset (bottom left): \V\ as a func- 
tion of substrate thickness for E s — 10 Pa. Other parameters: 
B c — 1 kPa, h c = 0.1 [im, l a — 6.3 ^m, L — 100 /im, u s = 0.4. 



within a length £ p from the edge of the cell layer. The 
penetration length £ p can be written as £ p = ^+1^, 



with = 



2B c h c h s 



the square root of the ratio of the 



cell's elastic energy to the elastic energy density of the 
substrate. This form highlights the interplay of focal ad- 
hesion stiffness and substrate stiffness in controlling spa- 
tial variation of stresses in the lateral (x) direction. The 
two act like springs in series, where the weaker spring con- 
trols the response. If Y a « 2 (i+v-s)h > ^ nen — ^ a anc ^ 
the stiff substrate has no effect. Conversely, if the focal 
adhesions are stiffer than the substrate, then l p ~ i s . For 
an elastic substrate with h s = 10 /im, v s = 0.4 and E s in 
the range 0.01 — 100 kPa, £ s lies in the range 0.2 — 17 /im. 
This leads to typical values of l p in the range 6.3—18 /xm 
for a cell layer of length 100 /im, consistent with exper- 
imentally observed traction penetration lengths on thin 
continuous substrates [10] and on micropillar posts [9]. 

Infinitely thick substrate. If h s » L, the substrate 
Green's function can be approximated as that of an 
elastic half plane, G(x) = — n £ E [7 + log (|x|/L)], with 
7 the Euler constant [3Tj. The solution of Eq. ^ 
with boundary conditions u(0) = cr(L) = can be ob- 
tained by expanding a(x) in a Fourier sine series as, 
a ( x ) — S^Li °V» sin (nnx/L) and solving the coupled al- 
gebraic equations for the Fourier amplitudes a n given in 
the Supplementary Material [57] . The effect of the nonlo- 
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cal elasticity of the substrate is controlled by yet another 
length scale £ soo = \J AB ^e c ^ that can be obtained from 
the length l s introduced in the case of thin substrate by 
the replacement h s — > L and (1 + v s ) — > 2/ir. This high- 
lights the known fact that cells or cell layers only "feel" 
the substrate up to a thickness comparable to their lat- 
eral size L. For parameter values quoted in the preceding 
paragraphs, £ soo takes values between 0.35 — 35 /zm for E s 
in the range 0.01 — 100 kPa, indicating that the thin/thick 
substrate crossover, although not observable in isolated 
cells, should be seen experimentally in cohesive cell layers 
where the lateral extent can exceed 100 /im. The cellu- 
lar stress and substrate displacement profiles obtained 
numerically by summing the Fourier series are shown in 
Fig. [3] (top) . The lateral variation of stresses is now con- 
trolled by the length scale £ p = £ 2 SOO . One con- 
sequence of nonlocal substrate elasticity is that the sub- 
strate deformation shown in the top right frame of Fig. [3] 
extends outside the region occupied by the cell layer, indi- 
cated by the two vertical dashed lines. The profile of the 
local traction stress displayed in Fig. [3] (bottom frame) 
shows that the traction stress is localized near the edge 
of the cell layer and its magnitude increases with sub- 
strate stiffness. The inset to Fig. [3] (bottom right) shows 
the magnitude of the net contractile moment defined as 
V = J_ oo dxxT(x). This quantity is negative, as ex- 
pected for contractile systems. Its magnitude increases 
with E s at a rate consistent with experiments, with a 
25% rise in \V\ upon increasing the substrate stiffness by 
40% [32], and saturates for very stiff substrates. 

Substrate of Finite Thickness. Finally, we consider a 
substrate of finite thickness, h s . The calculations are car- 
ried out using the approximate Green's function given in 
Eq. (J2J , with the replacement ( c q —> L. The variation 
of the net contractile moment with h s for E s = 10 Pa 
is shown in Fig. [3] (bottom left inset). As seen previ- 
ously in experiments [12] . \V\ drops sharply with increas- 
ing substrate thickness, quickly reaching the asymptotic 
value corresponding to infinitely thick substrates. Thin- 
ner substrates are effectively stiffer than thick ones, in- 
ducing larger contractile moments. Our analysis suggests 
a general expression for the penetration length £ p that in- 
terpolates between the thin and thick substrates limits, 

. / B c h c B c h c , 

^ = V^ + ^e ff . (6) 

Stress penetration is controlled by a substrate layer of 
effective thickness h„ = , „ A , — r + — given by the 

crt /i s 27r(l+z/ s ) L G J 

geometric mean of the actual substrate thickness h s and 
the lateral dimension L of the cell or cell layer. If h s << 
L, then h c s ~ 2?:h s (l + v s ) and stress penetration is not 
affected by cell layer size, as in the experiments of [TU] . 
On the other hand, if h s >> L, then cells only feel the 
effect of the substrate down to an effective depth L. 



Discussion. In summary, we have examined the de- 
pendence of traction stresses in adhering cell layers on the 
mechanical and geometrical properties of the substrate. 
Using a generic non-local model, we provide analytical 
results for the effect of cell and substrate properties on 
the stress penetration length, that can be tested in ex- 
periments. Although the analysis presented here is re- 
stricted to one dimensional layers, isotropic planar cell 
layers with spherical symmetry can also be considered 
analytically [33] . with similar predictions for the depen- 
dence of traction fields and their moments on substrate 
mechanical and geometrical properties. The scaling of 
traction moments on cell layer size is, however, different 
in two dimensions [10]. The model can be extended to 
incorporate the effects of cell polarization, spatial varia- 
tions in contractility, heterogeneities in the cell layer or 
anisotropic elasticity of the substrate. 
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